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Abstract 

A one-dimensional tight-binding Hamiltonian describes the evolution of a 
single impurity interacting locally with N electrons. The impurity spectral 
function has a power-law singularity A(oj) oc| uj — u>q with the same ex- 

ponent (5 that characterizes the logarithmic decay of the quasiparticle weight 
Z with the number of electrons N, Z oc N~@. The exponent (3 is computed 
by (1) perturbation theory in the interaction strength and (2) numerical eval- 
uations with exact results for small systems and variational results for larger 
systems. A nonanalytical behavior of (3 is observed in the limit of infinite 
impurity mass. For large interaction strength, the exponent depends strongly 
on the mass of the impurity in contrast to the perturbative result. 
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I. INTRODUCTION 



Anderson studied the effect of a static impurity potential on conduction electrons in 
metals and showed that the groundstate of the electrons is strongly renormalized by the 
local potential of the impurity and has an overlap with the unperturbed state, or quasipar- 
ticle weight Z, vanishing as Z oc with increasing number of electrons N. This effect, 
known as the Anderson orthogonality catastrophe, has its origin in an infrared singularity 
due to shake-up processes of the electron sea in the presence of the impurity potential |2| 
and signals the failure of the quasiparticle picture to describe the low-energy excitations. 

The infrared singularity also affects the optical properties of metals. The core-level hole 
created by an X-ray disturbs the conduction electrons similarly to an impurity potential. 
On one hand, the X-ray photoemission spectrum is asymmetrically broadened above the 
threshold @^]. On the other hand, the X-ray absorption spectrum has a strongly enhanced 
threshold, the so-called Fermi-edge singularity 0. 

These singularities in the optical spectra apply for a static core hole, i.e., an infinite- mass 
hole (or impurity). For a finite-mass hole and an isotropic band dispersion, the infrared sin- 
gularity does not occur in three dimensions because the hole recoil strongly restricts the 
number of low-energy excitations |||7| ; as a consequence, the edge singularities disappear. 
In one dimension, however, the infrared singularity persists even for a finite- mass hole ||. 
The observation of an enhanced threshold in UV-absorption spectra of doped semiconduc- 
tor quantum wires was interpreted as a Fermi-edge singularity || and stimulated renewed 



interest in the one-dimensional problem [10|JTl 



The present work studies the Anderson orthogonality for a finite-mass impurity in a 
simple one-dimensional model and focusses on the impurity recoil. Although the infrared 
singularity occurs for both an infinite- and a finite-mass impurity in one dimension, the recoil 
plays an important role on the critical exponent f3 which does not extrapolate to the static- 
impurity value in the infinite-mass limit. The dependence of the exponent on the impurity 
mass is investigated analytically by perturbation theory and numerically by a variational 
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approach in the nonperturbative regime. The results are eventually compared to analytical 
calculations from Ref. [O] which studied the heavy-mass and strong-coupling regime using 



a path-integral formalism (13]. 



The next section presents the model and the results. In section a perturbation 
analysis of the impurity spectral function is performed along the line of Ref. || . Section [TV] 
calculates numerically the critical exponent using a variational approach. 

II. MODEL AND RESULTS 

This section presents the model and summarizes the known results on the impurity 
spectral function and the Anderson orthogonality catastrophe. At the end of the section, 
the main results of this work are briefly described. 

The model describes a single impurity and N spinless electrons moving on a chain of 
L sites with periodic boundary conditions. Within a tight-binding approximation with 
nearest-neighbor hopping, the band energies are — 2th cos k and — 2t cos k for the impurity 
and the electrons, respectively. Further, the impurity and the electrons feel an interaction 
U when they sit on the same site. Although this study is restricted to a single impurity, it 
is convenient to write the Hamiltonian H in second-quantized form with creation operators 
cf for an electron on site i, and d\ for an impurity: 

L L L 

H = -tJ2 (cta+i + h - c ) - t h J2( d tdi+i + h.c.) + Uj24d4ci. (1) 

i=l i=l i=l 

The interaction is attractive in order to describe a hole in a valence band. For this particular 
model, however, the repulsive and attractive cases are related by a particle-hole transforma- 
tion for the electrons Cj = (— l) J 'cj. In the rest of the paper all the results are presented for 
the repulsive case; the corresponding results for the attractive interaction are obtained by 
the transformation p — > 1 — p where p = N/ L is the density of conduction electrons. 

The spectral function of the hole (or impurity) A(q,u) describes the photoemission re- 
sponse within the sudden approximation, i.e., neglecting interaction with the outgoing elec- 



tron ||14|| . It has a spectral decomposition in terms of eigenstates, | ip n ), and eigenenergies, 



E n , of the Hamiltonian (|T]) in presence of the impurity, and of the groundstate wavefunction, 
\ 4>o), and energy, E , in absence of impurity: 

A(g,u) = -1mG q (u) = I I 4 I I' 5 ( u ~ E n + Eo). (2) 

7T — * 

In the spectral decomposition of A(q = 0,u), the spectral weight of the groundstate Z =\ 
(if)o | dg =Q | <po) | 2 may remain finite in the thermodynamic limit, giving rise to a quasiparticle 
peak in the spectral function. This is the usual situation when the quasiparticle picture 
applies. 

Static impurity (th = 0). The interaction causes the quasiparticle picture to break down 
IJ: the spectral weight scales to zero with increasing number of fermions iV as Z oc iV _/3 , the 
groundstate of the interacting system being orthogonal to the quasiparticle state & \ 0o) f° r 
iV — > oo. This is known as the Anderson orthogonality catastrophe. The exponent is related 
to the phaseshift 6f of an electron at the Fermi energy scattered off the static impurity: 

P(t h = 0) = (5 F /7t) 2 . (3) 

For the present model, the phaseshift depends on U and on the density of states at the Fermi 
energy N F = l/(27rsin /cp): 5p = — arctan {ixUNp). Furthermore, the spectral function has, 
instead of a quasiparticle peak, a power-law singularity at the threshold ujq = Eq — Eq with 
the critical exponent 1 — j3(t h = 0) [|J: 

« I (4) 

| UJ — UJq 1 1 p 

This singularity in A{oS) is observed in X-ray photoemission of metals [[|. 

Finite mass impurity (th > 0). While an infinite-mass impurity acts as an external 
potential on the electrons, the impurity recoil further complicates the many-body problem. 
Despite this complexity, the eigenenergies and eigenstates of H are known exactly for the 
special case th — t |15|]. Using the exact solution, the spectral function for q ~ is computed 
in Ref. jnj. It has no quasiparticle peak because of the Anderson orthogonality catastrophe 
and has a power-law singularity with an exponent: 
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(3(t h = t) = 2 (5' F /n) 2 (5) 

The exponent is given by the phaseshift of a single electron at kp scattering off a finite- 
mass impurity 5' F = — arctan (irUNp /2). Notice the similarity in the exponents for th — t 
and th = 0, which are both expressed in terms of phaseshifts. The phaseshifts, however, 
differ since the number of states contributing to the Anderson orthogonality is reduced from 
UNp to UNp/2 between a static and finite mass impurity, respectively. The origin of this 



difference is discussed in Section |T| for the perturbative results. 

The present work presents calculations of the exponent P{th) for different hopping param- 
eters th, interaction strengths U, and electron densities p. The exponent (3{th) is computed 
(1) analytically using perturbation theory in U and (2) numerically in the nonperturbative 
regime using a variational wavefunction proposed by Edwards ||17|| . The exponent is ex- 
tracted numerically from finite size results using a precise scaling law for Z as a function of 
N and a numerical fit of the data as N — ► oo. The main results of this study are summarized 



now. Sections |TJ and [TV] will give a detailed description of the perturbation calculations 
and the numerical simulations, respectively. 

Mass dependence of the exponent. The perturbative results indicate that the exponent 
P, to order (U/t) 2 , is independent of th for th > and equals (UNp) 2 /2. For a finite U, 
however, the exponent P(th) calculated numerically does depend on t h and its dependence 
increases with increasing U, as illustrated in Fig. [I] that shows the exponent normalized 
to its value for th = t, P{th)/P{th — t), as a function of th- In the strong coupling limit, 
U = oo, (3 varies quasi linearly with th- 

Discontinuous exponent in the heavy-mass limit. The perturbative calculations predict 
a discontinuous exponent in the limit of a flat impurity dispersion: lim^^o Pith) = Pith = 
0)/2. This discontinuity is due to the irrelevance of backscattering processes for th > 
because of the finite recoil energy involved. The numerical results in Fig. p] illustrate the 
nonanalyticity for a finite U: when th — > the exponent does not extrapolate to the static 
value which is indicated by the filled symbols at th = 0. The numerical results can be 
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compared to calculations of A. Rosch and T. Kopp |TI| for the exponent in the heavy-mass 



and strong-coupling regime. Their analysis, based on an effective action for the long-time 
behavior of the impurity propagator, also predicts a discontinuity lim th ->o (3 (th)/ P(t) = a ^ 
(3(0) /(3(t) and gives at half filling a = 1/4 while at third filling a = 19/56 ~ 0.339. The 
numerical results in Fig. [IJ suggest a ~ 0.25 for both half and third filling. The result 
at half filling is therefore in good agreement with their prediction while at third filling the 
value of a is significantly smaller. Note however that the present work relies on a variational 
approach. 

Crossover behavior in the heavy-mass limit. A detailed analysis of the perturbative 
results for th <C t reveals a crossover in the scaling behavior of In Z as a function of In N, 
as shown in Fig. [|. While \nZ closely follows the static-impurity behavior with a slope 
f3{th = 0) for a number of fermions smaller than a crossover value N c , it adjusts to the 
finite-mass behavior only for iV > N c . Further N c diverges as t/th for th — > 0. Therefore, 
the discontinuity is only an asymptotic result, valid in the limit N — > oo. 

III. PERTURBATION THEORY 

This section evaluates the impurity spectral function and spectral weight Z perturba- 
tively in U and for an arbitrary hopping parameter th in one dimension, following Ref. || 
that computes A(q = 0, to) for the equal-masses case (th = t) in connection with the stability 
of the ferromagnetic state in the Hubbard model. 

Spectral weight. In the perturbative expansion, the first terms that renormalize the 
ground-state wavefunction correspond to the creation of a single particle-hole pair within 
the Fermi sea by an impurity of momentum q — 0. The excitation energy is: Ae(ki, k^) = 
2t cos ki — 2t cos k 2 + 2th — 2th cos(fci — k 2 ). The spectral weight has a cumulant expansion 
0] that, up to second order in U, involves only these excitations: 

, „ _ fU\ 2 v 9(fc F -|fa|)9(N-M 
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The sum over ki,k 2 diverges logarithmically with increasing number of electrons N. As 
shown in the Appendix, a large A" expansion gives : 

\nZ = -{3(t h )\nN + a (t h ) + ai (t h )/N + 0(l/N 2 ). (7) 

The logarithmic term dominates for large A" and gives rise to the Anderson orthogonality. 



The finite-size corrections are used in Section IV for the numerical study of (3. 



The main result of this section is the evaluation of (3(t 



h 



13(0) = (UN F ) 2 for t h = 

f3{t h ) = \{UN F ) 2 forf„>0 (8) 

The exponent (3{t h ) is independent of t h as far as t h > 0. The infrared singularity is caused 
by forward-scattering processes with small momentum transfer \k 2 — k\\ <C kp, which are 
gapless excitations for any th- The hopping th is irrelevant since the impurity recoil energy 
2t/ l (cos(/ci — k 2 ) — 1) is negligible as compared to the particle-hole energy 2t(cos k 2 — cos ki). 

Further, the exponent has a discontinuity in the heavy- mass limit lim^^o fl(th) — P{th = 
0)/2. The difference between the infinite- and the finite-mass exponents is simply related to 
scattering of one electron from one side of the Fermi surface to the other. These so-called 
backscattering processes, which involve a large momentum transfer \k 2 — k\\ ~ 2k f, do 
not contribute to the infrared divergence for th > since the impurity recoil opens a gap 
2th{cos2kF — 1). For th — however, both backscattering and forward-scattering processes 
are gapless. The number of low-energy excitations contributing to the infrared singularity 
is thus reduced by a factor of 2 for a finite-mass impurity as compared to its value for the 
static impurity. 

The discontinuity in the exponent is an asymptotic result valid only for — > oo. For 
a finite number of fermions and a large but finite mass < th t however, InZ has a 
crossover as a function of In N, illustrated in Fig. |2] where the spectral weight is computed 
numerically from (^|). The logarithm of the spectral weight has the slope (3(th) only for 
a number of electrons larger than a crossover value N c , while for a small it follows the 



static- impurity behavior with a slope (3{th = 0). The dashed lines indicate the asymptotic 
behaviors — InZ = (3(th) InN — ao(th), with the analytical (3 from (|8]) and a® fitted to the 
value of In Z for the largest size. The intercept of the asymptotes gives the crossover size 
N c , which agrees very well with the estimate iV c ~ 0.3244t/th, presented in the Appendix. 
Notice that N c diverges as t/th for t h — > and the asymptotic regime is reached for a larger 
number of electrons the lower th- 

Spectral function. The spectral function is computed only at q = where it has a power 
law singularity. The propagator is: 

G q=0 (r) = - i (0o | d q=0 exp(-tHr)dl =0 | O > e l ^ T Q{r). (9) 

The propagator has also a cumulant expansion and is written in terms of the density of 
particle-hole excitations S(u) and a renormalized impurity energy eg: 

G q=0 (r) = - i exp 

S(oj) = — ^ T f e(|fci| - k f )e(k f - \k 2 \)S(uj - Ae(h, k 2 ))dhdk 2 . 
(2n) J -it J--* 

The density S has a linear frequency dependence for small u: S(u) = (/3{th)/U 2 )u where 
/3(th) is the exponent of the Anderson orthogonality in (^|). This linear behavior determines 
the low-frequency spectral function which has a power-law singularity at threshold with the 
exponent 1 — (3 0] : 

a/ r\ \ sin(7T/3)r(l — (3) , 

TT{LU-e ) L P 

This is however only an asymptotic result for frequencies smaller than a cutoff W. For 
th = 0, the cutoff is of the order of the Fermi energy: W ~ 2t(l — cosfcp). For th > 
however, the linear behavior of S holds only for frequencies smaller than the impurity recoil 
energy and the cutoff is given by W ~ min(2t(l — cos kp), 2^(1 — cos2fcp)). For a heavy 
impurity the cutoff is of the order of the impurity recoil energy £^(1 — cos2/c^) rather than 
the Fermi energy and the asymptotic result is valid only in a very narrow frequency range. 
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S{u>) 



exp( 



-lUJT) 



du 



(10) 



Furthermore, the density of excitations exhibits a crossover similarly to the spectral weight. 
This might give rise to a crossover in the spectral function as well. 

In summary the exponent (3 in (|8[) characterizes the power-law singularity of the spectral 
function and the logarithmic decay of the quasiparticle weight. The exponent does not 
depend on the mass of the impurity except in the static limit = and it has a discontinuity 
at th = 0. Notice that the perturbative results agree with the the small-?/ expansions of the 
exponent for th = and t^ — t. 



IV. NUMERICAL STUDY 

This section presents a numerical study of the exponent (3{th) based on a variational 
approach. The variational predictions for the energy and correlation functions are com- 
pared to results from Lanczos exact diagonalizations and Projection Quantum Monte Carlo 
simulations. Then the variational calculations are used to extract the exponent. 

The variational wavefunction was originally proposed for the single spin-flip problem in 
the 2D-Hubbard model in reference to the stability of the ferromagnetic state |17j . It was also 
used to study numerically the quasiparticle weight in two dimensions ||18||. Furthermore, this 



variational approach is equivalent to the approximation used in Ref . JTTJ . In one dimension, 
the variational class of wavefunctions contains all the eigenstates of the model (|I|) for t% = t, 
as shown by Edwards |T7|]. The variational approach is thus expected to include much of 



the relevant correlations even for th ^ t. 

In the reference frame comoving with the impurity a wavefunction | ty q ) of total mo- 
mentum q is represented by a function f(ji, . . . ,Jn) depending only on the positions of the 
electrons: 

I *,) = iE^< E A/, ./n ./v-./n)-...- x 0). (12) 

v io=i ii,-Jiv=i 

The variational ansatz for / is a determinant of single-particle wavefunctions <p m : 

f(ji,---,3N) = 7/^fy det [0mC7'z)] m , i= i ,...,n- ( 13 ) 



The expectation value of the energy is: 



(H) = E£i [-t?:}^m)M + 1) + cc.) + C/ 1 0,(0) 
— £/, [exp(— zg) det(S) + cc] . 

Smn=E€(i + 1 M)' ( 14 ) 
3=0 

The variational parameters <fi m (j) are found by minimization of the energy using a steepest 
descent algorithm. If one chooses as starting 4>i{j) the exact solution for th = t, convergence 
is reached after a relatively small number of iterations even for hopping parameters very 
different from t. 

Comparison of variational results to Lanczos exact diagonalizations and Projection Quan- 
tum Monte Carlo simulations. Only small systems are accessible by exact diagonalization 
because the dimension of the Hilbert space increases very rapidly with the number of lattice 
sites. For the half-filled band, the biggest closed-shell system studied has L = 22 whereas 
for p = 1/4 it has L = 20. For larger systems, Monte Carlo simulations are required. 

Projection Quantum Monte Carlo gives a statistical estimate of ground-state expectation 
values of observables O by a projection of a trial wavefunction | ip T ) onto the groundstate 
with the operator exp(—r]H) |i~9| : 





exp(—i]H)0 exp(—i]H) 


V»r) 




| exp(-r]H)\il) T }\\ 2 



(ipo | O | = lim^oc . (15) 



The algorithm used here closely follows Ref. | 20| . The imaginary time evolution exp(-rjH) 



is performed sequentially for small time intervals At and a Trotter decomposition is used for 
the kinetic-energy and the interaction terms in the Hamiltonian. Eventually the two-body 
term is represented by discrete Hubbard-Stratonovitch fields that mediate the interaction. 
For the present calculations a rather large rj = 15 is necessary in order to converge the 
relevant correlation functions, and the results are extrapolated to At — ► using several 
values of the time interval. 

The energy is computed for the groundstate, which has zero total momentum. The 
relevant quantity is the correlation energy e c = E — E — Up where E and E are the 
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interacting and noninteracting energy, respectively. Notice that e c is of order 1 while Eq 
is of order N. Table | gives e c for several system sizes at half and quarter filling and for 
an interaction equal to the bandwidth U = At and a hopping parameter th = 0.51 The 
variational ansatz is not exact for th ^ t since the variational energy departs significantly 
from the exact energy. Nevertheless it remains always very close to the ground-state energy. 
The relative difference between the exact and the variational energy remains smaller than 
0.1% even when extrapolating to the infinite-system limit with a 1/N scaling law. For 
the sake of comparison the correlation energy for the unrestricted Hartree-Fock solution is 
— 0.4380i at quarter filling, a value much larger than the variational energy — 0.64131 

The k = component of the momentum distribution function of the impurity n(k = 0) 
has the same scaling law as Z when the Anderson orthogonality catastrophe occurs. Both 
Z and n(k = 0) are computed in order to test the relevance of the variational calculations 
to extract the exponent j3. Figure |3| shows the ratio of the exact momentum distribution 
n(k = 0), calculated by either exact diagonalization or Quantum Monte Carlo depending 
on the system size, to the variational estimate n(k = 0) var as a function of 1/N at quarter 
and half filling, for U = At and t^ = 0.51 If the variational result reproduced correctly the 
scaling behavior of n(k = 0) with N, the ratio n(k = 0)/n(k = 0) mr should be constant 
as a function of 1/N. At half filling, the ratio remains indeed always close to 1 and the 
variational approach seems to correctly describe the scaling behavior. At quarter filling 
however, the comparison relies mostly on Monte Carlo simulations where the estimate of 
n(k = 0) has strong statistical fluctuations for large system sizes, large autocorrelation times 
having an essential contribution to the error bars. The comparison to the variational results 
is therefore delicate but still the ratio remains close to 1 within the error bars. Finally 
exact-diagonalization results for Z are presented in Fig. ^for U — oo, p — 1/2 and different 
th- The slope of In Z as a function of IniV is an estimate of the exponent (3 and both 
exact-diagonalization results and variational calculations are in good agreement. These last 
results illustrate the robustness of the variational wavefunction even in the strong-coupling 
regime. 
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Finite-size scaling analysis of the exponent. The finite size corrections to the Anderson 
orthogonality catastrophe in (|7|), which were derived from the perturbative analysis, are used 
now to extract numerically the exponent. It should be stressed that these 1/N corrections 
contrast with the slowly decaying lnlniV corrections expected for the paramagnetic phase 
of the 1-D Hubbard model J2T] . 

The exponent is extracted from numerical data on finite systems by scaling the slope 
of In Z as a function of IniV: the spectral weights Z\ and Z2 are computed for systems 
with a number of fermions N\ and N2, respectively, at a fixed density and a fixed ratio r = 
N1/N2] from the perturbative analysis in (0), the slopes (In Z\ — hi Z%)/ (In N2 — In N%) have 
a polynomial expansion in 1/N, N being the mean number of electrons N = (JVi + iV 2 )/2: 



The exponent (3 is estimated by a numerical fit of the slopes as N —>■ 00. In practice, it is 
not possible to keep the ratio Ni/N 2 exactly fixed while increasing N but the number of 
electrons can be adjusted such that Ni/N 2 approaches a fixed value for large N. All the 
results presented here have been obtained for Ni/N 2 — 0.7. 

The scaling procedure is tested both in the perturbative regime for different th and at 
finite U for th = t, where f3 is known analytically. Figure [5] illustrates the scaling procedure 
for the perturbative regime where InZ is computed from @. The inset shows that the 
scaling procedure is well behaved even for very large N. Furthermore, the scaling analysis 
is essential to determine the dependence of (3 on th since the finite-size corrections differ in 
both sign and magnitude for different th- Note also that the non-linear terms in 1/N become 
increasingly important with decreasing th- The spectral weight Z is computed numerically 
for th = t and different U using the Bethe's ansatz wavefunction (|T3| ) and the results are 
fitted with a third-order polynomial in 1/N. The relative accuracy of the fitting procedure in 
extracting j3 remains of the order of 10 -4 even for a strong interaction U = 8t, as illustrated 
in Fig. ^. Although the scaling behavior was derived in the perturbative regime, it seems 
to hold for any interaction strength U. 




(16) 
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For a finite U and th 7^ t, the exponent is not known and its evaluation relies on the 
variational approach and the finite-size scaling analysis. Figure |7| illustrates the scaling 
procedure in the strong-coupling regime U = 00, for p = 1/2 and different th- For th > O.lt, 
the exponent, which is obtained by an extrapolation of the data as iV — ► 00, depends only 
slightly on the order of the polynomial used in the fit. For th = O.lt, however, the fitting 
procedure is not well-behaved and a relative error of a few percents is expected in the 
extraction of (3. A precise investigation of the heavy-mass regime th <S t would require the 
study of even larger systems due to the important nonlinear corrections in 1/JV. 

Mass dependence of the exponent. The study of Pith) is based on the variational approach 
and the scaling analysis presented above. Figure [1] presents the exponent normalized to its 
value at t h = t, /3(th)//3(t h — t), as a function of t h /t. For a finite U, the exponent does 
depend on the hopping parameter t h , in contrast to the perturbative result Pith) /Pith = 
t) = 1 indicated by the dotted line. Furthermore, this dependence increases with increasing 
U and in the strong-coupling regime U = 00 the exponent varies quasi-linearly with th/t. 
For small th and at half filling, the exponent only slightly departs from the linear behavior. 
Furthermore, the exponent depends on the density since the data for U = 00 and p = 1/2 
significantly differ from the exponents for p = 1/3. This contrasts with the exact exponent 
for th = and th = t which are independent of the density. However, the limiting value 
lim th ^o Pith) seems independent of the density. 

The numerical results demonstrate the discontinuity of the exponent at th = for a finite 
U. The numerical data, indeed, do not extrapolate to the exact results for th = indicated 
by the filled symbols in Fig. |]. A precise extraction of the limiting value lim^^o Pith), 
however, would require the investigation of larger systems. 

The occurrence of the nonanalyticity at th = is not surprising since the translational 
symmetry is broken at this point. Still the role played by the recoil of the impurity in the 
discontinuity of the exponent is not clear. The perturbation calculations indicate that the 
discontinuity at th — is due to the irrelevance of the backscattering processes whenever 
th > 0; yet a simple argument can persuade us that this is true only in the small U regime. 
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Indeed let us assume first backscattering to be responsible for the nonanalyticity for all U. 
One can devise an effective model for the heavy-impurity limit (t^ <C t) where the impurity 
is considered static but the interaction with the electrons is restricted to forward scattering : 

H = -2tJ2 cos k c{c k + Y c k c k'- (17) 

k kk'>0 

In this picture, the only effect of the impurity recoil is the restriction to forward-scattering 
processes. Since the potential is static, the exponent is expressed in terms of phaseshifts. A 
calculation of the phaseshift gives the same exponent as the result for th — t in ([5]). This is 
not compatible with the numerical results for <C t. Therefore forward-scattering processes 
alone cannot account for the discontinuity and backscattering has to be invoked. 

V. CONCLUSIONS 

This study of the one-dimensional Anderson orthogonality catastrophe, combining ana- 
lytical and numerical calculations, has focussed on the effect of the impurity recoil. 

The numerical study requires a finite-size scaling analysis since the Anderson orthogonal- 
ity catastrophe results from a logarithmic decay of the quasiparticle weight with the number 
of fermions. The present work shows, however, that a reliable numerical analysis of the 
quasiparticle renormalization can be achieved if a precise scaling hypothesis is established 
and large enough systems are accessed. 

Within perturbation theory, the infrared singularity that signals the orthogonality catas- 
trophe occurs for an impurity band of any dispersion. Still, there is a discontinuity between 
the zero-bandwidth exponent and the finite-bandwidth exponent. This discontinuity is re- 
lated to the impurity recoil which opens a gap in the spectrum of particle-hole excitations 
for backscattering of one electron from the Fermi momentum kp to —kp. 

Outside of the perturbative regime, the numerical analysis demonstrates the discontinu- 
ous behavior of the exponent and agrees at half filling with a study of the heavy-mass and 
strong-coupling regime in Ref. J12|| . At third filling however, the numerical results differ 
from the analytical prediction. 
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The discontinuity of the exponent is an asymptotic result valid only in the limit of an 
infinite system. For a finite number of electrons and a heavy but finite-mass impurity the 
quasiparticle weight has the same logarithmic behavior as a static impurity up to a critical 
number of electrons N c where the former weight shows a crossover to the true asymptotic 
decay for a finite- mass impurity. Furthermore, N c diverges with increasing mass of the 
impurity. This crossover is also expected in the low-frequency behavior of the spectral 
function. 
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This appendix derives the scaling law fl7|) from the perturbative expression of the spectral 
weight Z and estimates the crossover size N c . The calculations are presented for the half- 
filled band in details and the density dependence of N c is briefly discussed at the end. 

Forward scattering processes. An electron with momentum k\ is scattered into an empty 
state with momentum k 2 such that kik 2 > 0. The infrared singularity is caused by excitations 
around the Fermi momentum kp = 7r/2 whose energy vanishes linearly with the momentum 
transfer ki — k 2 . For ki and k 2 > 0, the momenta are written as k\ = kp — 2 / im 1 /L, 
k 2 = kp + 27m 2 /L and the excitation energy is: 



ACKNOWLEDGMENTS 



APPENDIX: 




(Al) 



15 



While the linearization of the energies allows an exact calculation of the exponent /3(th), 
it provides only an estimation of the remaining terms in (|7|). The contribution of forward- 
scattering processes to the spectral weight is: 

U\ 1 1 N/2 N/2 1 

TJ jf >0 Ae(k u k 2 ) 2 = '2 {UNF)2 ^ ^ 1 (n 1 +n 2 y(l + nt h (n 1 +n 2 )/(2Nt)r 

(A2) 

The fraction in the sum is expanded in four different terms: 

1 1 1 

+ 



(n l + n 2 ) 2 (l + ^t h {n 1 + n 2 )/{2Nt)) 2 (m + n 2 ) 2 (m + n 2 + 2Nt/(>nt h )) 2 

+TT1 — „ u u (A3) 



Nt(rn + n 2 ) Nt(n 1 +n 2 + 2Nt/(7rt h )) 



The first term gives the logarithmic divergence |22 

if = m W + 1 + - 2.n2 + I + (i) . (A4, 

where C ~ 0.5772 is Euler's constant. The contributions of all other terms remain finite for 
N — > oo. For the second and fourth terms in( |A3| ), the discrete sum can be replaced by an 
integral. Further the finite-size corrections are of order 1/N. The extraction of the finite-size 
corrections for the third term, however, requires the evaluation of the discrete sum: WM: 



1 N/2 N/2 1 / 1\ 1 / 1 \ 

-Y T =lB2+(ln2-i)- + 0(-= 7 ). (A5) 

Backs cattering processes. The initial and final momenta have a different sign, k\k 2 < 
and, for k\ > and k 2 < 0, the momenta are written as k\ = kp — 2^n\/L and k 2 = 
—kp — 27m 2 /L. The linearized excitation energy is: 

M^ 2 ,)=4(4 + ^!±^). (A6) 

For th — 0, backscattering processes are gap less excitations and have the same contribution 
to Z as forward-scattering processes, as given in ( |A4j ). For th > 0, however, they have a 
gap and their contribution, which is finite, is computed by an integral representation of the 



sum. 
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Finally, all the terms, for both forward-scattering and backscattering, give finite size 
corrections of order 1/N. Furthermore the coefficient a® in (|7|) is evaluated explicitely in 
order to extract the crossover size iV c : 



«o(th = 0) = -/3(0)(l + C-21n2). 



(A7) 



a (t h > 0) = -P(t h ) 



1 + C + In 



2t(irt h + 2t) 7it h nt h + 2t (vrt + 4t 



+ ^ln 



+ In 



(nt h + 4t) 2 t irt h + At 8t h (irt + 2t h ) 
The crossover size iV c is the solution of the equation: 

- 0(0) lnN c + oto(O) = -P(t h ) lniV c + a (4) 



• (A8) 



(A9) 



The crossover occurs for ^ -C t where a (t h ) is diverging logarithmically. As a consequence 
iV c is proportional to th/t: 



nt t 
N c = — exp(-l - C) ~ 0.3244—. 
2th £/i 



(A10) 



Arbitrary density. The scaling law (0) applies to any density. Furthemore the crossover 
size is proportional to t/th and diverges as 1 — p in the limit 
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FIG. 1. Exponent (3(th)/ P{th = t) normalized to the analytical value (3(th = t) in (||) as a 

function of the mass ratio th/t for different interaction strengths U and densities p. The filled 

symbols at th = are the exact results in(|3|) for the static impurity. The open symbols for 

t/j > are the exponents extracted as ./V — > oo from the variational results for a sequence of finite 

systems with N1/N2 = 0.7, the largest size being Ni = 159 and N\ = 121 at half and third filling, 

respectively. 
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FIG. 2. Logarithmic decay of the quasiparticle weight Z with number of electrons N from the 
perturbative result in (|6|) at half filling and for different impurity hopping parameters t^. The 
dashed lines are the asymptotic behaviors with the analytical exponent f3(th) in The arrows 
indicate the estimate of the crossover size N c from ( |A10 ). 
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FIG. 3. Ratio of the momentum distribution function n(k) at k = 0, calculated with Lanczos 
diagonalization or Projection Quantum Monte Carlo, to its variational estimate n(k = 0) var as a 
function of 1/N for half and quarter filling and for U = At, th = 0.5i. The circles are the Lanczos 
results and the squares the Monte Carlo results with the corresponding error bars. The left axis 
refers to half filling and the right axis to quarter filling. 
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FIG. 4. Comparison of exact spectral weight Z from Lanczos diagonalization to its variational 
estimate as a function of number of electrons N in log-log plot forL r = oo,p = l/2 and different 
th- The open symbols are the variational results and the filled symbols Lanczos results. 
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FIG. 5. Slopes dln(Z)/dln(N) = (\nZi — In Z2) / (In N\ — \1\N2) normalized to the analytical 
value of f3(th) in (|8|) as a function of the inverse mean number of electrons 2/{N\ + N2) in the 
perturbative regime and for different hopping parameters t^. The spectral weight Z is computed 
numerically from (^). The dashed lines are numerical fits of the data as polynomials in 2/(iVi +./V2). 
The inset is a blow-up around the origin. 
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FIG. 6. Slopes dln(Z)/d ln(iV) normalized to the analytical value of f3(th = t) in @ as a 
function of the inverse mean number of electrons 2/{N\ + N2) for th = t and different interaction 
strengths £/ = 2,4, 8i at half filling and for N1/N2 = 0.7. The spectral weight Z is computed 
with the Bethe's ansatz wavefunction in The dashed lines are numerical fits of the data as a 
third-order polynomial in 2/(N\ + iV^). The inset is a blow-up around the origin. 



22 



N 

s 



• t h =0.5t 


U=oo p=l/2 


■ t h =0.25t 

" S7 U 
3 order fit 


Nj/N^O.7 


4 order fit 




* 












89 • 

61 


47 / 7 ' 'J 



0.000 0.010 0.020 0.030 

1/(N,+N 2 ) 

FIG. 7. Slopes d\n(Z) / d\n(N) normalized to the fitted (3{th) as a function of the inverse mean 
number of electrons 2/{N\ + N 2 ) for U = 00, different hopping parameters t/j = 0.5,0.25 and 
O.lt, at half filling. The spectral weight is computed with the variational wavefunction. The 
solid and dashed lines are numerical fits of the data as third-order and fourth-order polynomial in 
2/(N 1 + N 2 ), respectively. 
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TABLES 



N(p=l/2) 


Lanczos 


e c 
PQMC 


variational 


N(p=l/4) 


Lanczos 


e c 
PQMC 


variational 


5 


-0.9028 




-0.9025 


5 


-0.6102 




-0.6101 


7 


-0.9230 




-0.9227 


7 




-0.619(4) 


-0.6190 


9 


-0.9340 




-0.9336 


9 




-0.624(3) 


-0.6239 


11 


-0.9408 




-0.9404 


11 




-0.627(2) 


-0.6271 


13 




-0.945(2) 


-0.9451 


15 




-0.632(2) 


-0.6309 


oo 


-0.9727 




-0.9720 


oo 




-0.642(4) 


-0.6413 



TABLE I. Comparison of correlation energy, e c = Eq — Eq — Up where Eq and Eq are the 
ground-state energies of the interacting and noninteracting system, respectively, computed by Lanc- 
zos exact diagonalizations, Projection Quantum Monte Carlo (PQMC) and Edwards' variational 
wavefunction, as a function of the number of fermions N at half and quarter filling (p = N/L = 1/2 
and 1/4, respectively), for th = 0M, U = 4t. The last line gives the energy extrapolated for the 
infinite system both from either Lanczos or PQMC results and the variational results by a linear 
fit in 1/JV. 
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